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Abstract 



We study the influence of an external driving field on the coherence properties 
of a qubit under the influence of bit-flip noise. In the presence of driving, two 
^ , paradigmatic cases are considered: (i) a field that results for a suitable choice of the 

\^ \ parameters in so-called coherent destruction of tunneling and (ii) one that commutes 

■ with the static qubit Hamiltonian. In each case, we give for high-frequency driving 

, a lower bound for the coherence time. This reveals the conditions under which the 

\ external fields can be used for coherence stabilization. 
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qh| 1 Introduction 

> ' 

^ ! The experimental realization of one-qubit gates in solid state setups [1,2] and 

\ two-qubit gates in ion traps [3,4] and Josephson junctions [5] has demonstrated 

that these systems provide remarkable coherence properties although the goal 
of 10~^ errors per gate operation [6] has not yet been achieved experimentally. 
Theoretical studies of decoherence of two-level systems [7] have been extended 
to gate operations in the presence of an environment in Refs. [8,9, 10, 11, 12]. 
The unavoidable coupling to external degrees of freedom and the thereby 
caused decoherence still presents a main obstacle for the realization of a quan- 
tum computer. Thus, proposals for the stabilization of the coherence of qubits 
are particularly welcome. 

A variety of suggestions in this direction relies on the control of coherence 
by the influence of external fields. In particular, it has been proposed to use 
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the physics of the so-called coherent destruction of tunneling (CDT) for that 
purpose. CDT has originally been discovered in the context of tunneling in 
a driven bistable potential [13,14]. There, it has been found that a particle 
which is initially in the, say, left well of a symmetric bistable potential, can 
be prevented from tunneling by the purely coherent influence of an oscillat- 
ing driving field. This effect is stable against dissipation in the sense that the 
ac field also decelerates the dissipativc transitions from the left to the right 
well [15, 16, 17] thereby stabilizing coherence. For a preparation in a delocal- 
ized superposition, it has been found [10] that the coherent dynamics is not 
suppressed. 

A different proposal for coherence stabilization is the application of a sequence 
of TT-pulses that flip the sign of the qubit-bath coupling operator resulting in 
a so-called dynamical decoupling (DD) of qubit and bath [18,19]. A draw- 
back of this scheme is the fact that it eliminates only noise sources with a 
frequency below the repetition rates of the pulses. This clearly causes prac- 
tical limitations. However, these limitations may be circumvented by using a 
related scheme based on continuous-wave driving, i.e. one with a harmonic 
time-dependence, which allows higher driving frequencies. 

In this paper we investigate the coherence properties of a qubit coupled to 
an ohmic environment which we model by a spin-boson Hamiltonian [20]. 
We consider two different types of harmonic driving: First, one that leads 
to coherent destruction of tunneling and, second, one that corresponds to a 
continuous wave version of the dynamical decoupling pulses. In Section 2 we 
review the spin-boson model for weak dissipation and derive in Section 3 a 
Markovian formalism for driven dissipative qubits in the form of a basis-free 
master equation. The global decoherence will be qualifled by an upper limit 
for the decoherence rate that already gives a reliable estimate for the actual 
value. Subsequently, in Sections 4 and 5, we study the coherence control by 
oscillating external fields. 



2 Qubit with bit-flip noise 



The entire system of the qubit and the environmental degrees of freedom is 
described by the microscopic Hamiltonian 

H — i^qb + -f^coupl + -f^bath, (1) 

where the qubit is formed by a two-state system with a level splitting hA, 
thus 



2 



Below, we will in addition consider an oscillating driving field acting on the 
qubit. The environment is modeled by a bath of harmonic oscillators which 
couple linearly to the qubit, 

-f^bath = XI ^ialoi- (3) 

i 

-f^coupi = \(yx X ^Cj(a| + ai) , (4) 
^ i 

where Ui are the oscillator frequencies and hci the qubit-bath coupling energies. 
The bath couples to the qubit operator Ux thereby inducing bit fiips, i.e. 
incoherent transitions between the ground state and the excited state of the 
qubit. 

We complete the model by choosing an initial condition ptot(^o), of the Feynman- 
Vernon type, i.e., we assume that at time t — to the bath is in thermal equi- 
librium and uncorrelated with the qubit, 

Ptot(^o) = p{to) ® -Rbath,eq, (5) 

where p is the reduced density operator of the qubit and -Rbath.cq exp(— /Sifbath) 
is the canonical ensemble of the bath at the inverse temperature /? = 1/kBT. 



2.1 Dissipative qubit dynamics 

Prom the Liouville-von Neumann equation ihptot — [H, ptot] for the total den- 
sity operator one obtains for the reduced density operator of the qubit by 
standard techniques, the Markovian weak-coupling master equation 

i 1 /'°° / 

h AJo \ (6) 

+A{t) [a^, {a^it - T, t),p}]) , 

where {A, B} = AB + BA stands for the anticommutator. The notation 
X{t,t') is a shorthand for W{t,t')XU{t,t') with U being the propagator for 
the coherent qubit dynamics. Note that for later use, we have written the 
master equation already in a form which is also valid in the case of an explicit 
time-dependence of the qubit Hamiltonian. The dissipative equation of motion 
(6) depends on the bath through the symmetric and anti-symmetric correla- 
tion functions of the bath operator B{t) = J2i Q exp^iujit) + ai exp{—iuit)), 
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given by 



S(t) = I ({B(t),B(0)}) = - r dujJiuj) coth(^/3/2) cos(cur), (7) 

2 TT Jo 

A{t) = \ ([B(t),B(0)]) = -- da;J(a;)sin(a;T), (8) 

2 TT JO 

respectively. The angular brackets (■ ■ ■ ) denote the average with respect to 
the thermal equilibrium of the bath. We have introduced for the qubit-bath 
coupling the spectral density 

J{uj)^T:Y.cl5{u;-uJi). (9) 

i 

If the bath modes are dense, J{uj) becomes a smooth function. Within the 
present work, we will consider the case of Ohmic dissipation where J{u}) ~ 
I'KaLO exp(— ci;/a;c) and a is a dimensionless measure for the dissipation strength. 
The so-called cutoff frequency lOc is assumed to be the highest frequency of the 
bath. If the cutoff frequency is much larger than all relevant energy scales, the 
antisymmetric correlation function A{t) can be replaced by 27riQ;5'(r). Then, 
for the integral in the master equation (6), the part containing A(r) can be 
evaluated after integrating by parts to read 

— [a„{[i/V^,<7,],p}]. (10) 



To bring the master equation (6) into a more explicit form, we insert the 
Heisenberg operator 

^xi^ — Tjt) — CTx cos(Ar) + ay sin(Ar), (11) 

which is readily derived from its definition together with the qubit Hamiltonian 
(2). Performing the integration over r and uj and neglecting renormalization 
effects, which are small provided that aln{ujc/A) <S 1, yields for A <^ lOc, the 
Markovian master equation 

i r TTQ; A 

P= -^[Hqb,p] - -^[(^x, ['^x,p]]+i^^[<7x,Wy,p}]- (12) 

We will find below that the relaxation and decoherence processes are deter- 
mined by the rate [20] 

r = ls{A\ (13) 

where 

S{uj) = 27rauj coth j (14) 

is the power spectrum of the bath fiuctuations, i.e. the Fourier transformed of 
the symmetric bath correlation function (7). This master equation is a basis- 
free version of a Bloch-Redfield or Floquet-Markov master equation [21,22]. 
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Such an operator notation is preferable to a decomposition into the qubit's 
eigenbasis for ease of notation and, moreover, since it allows a more elegant 
computation of expectation values. The interpretation of (12) is that the first 
term of the right-hand side is responsible for the coherent dynamics, while the 
second and third term correspond to decoherence and relaxation. 

For the description of dynamics of a single qubit, it is convenient to map the 
density operator to the Bloch vector s = tr (cxp). It is straightforward to derive 
for the Bloch vector from the master equation (12) the inhomogeneous linear 
equation of motion 

k^-Ms + b, (15) 

where 



M 






-A 


A 


r 





r 






-vraA 



(16) 



In the weak dissipation limit, F <c; A, the matrix M has the eigenvalues 



1. 



-F±iA. 



(17) 



While the first eigenvalue describes the decay of the population of the excited 
state, the other two correspond to damped oscillations of the off-diagonal 
density matrix elements in the basis eigenbasis of the qubit Hamiltonian. This 
justifies for F the designation relaxation rate and the definition of the deco- 
herence rate F^ = F/2. For larger systems there is room for several relaxation 
(purely real eigenvalues) and decoherence rates (complex eigenvalues with non- 
vanishing imaginary parts). 



2.2 Entropy production 



Relaxation and decoherence are nonunitary processes, i.e. processes which 
transform a pure state into a mixed state. Since a quantum computer relies 
fundamentally on the coherence of the time evolution, such processes put an 
essential limitation and their influence has to be minimized. Thus, in order to 
follow the process of coherence loss, it is desirable to define proper coherence 
measures which preferably are independent of the chosen basis. A natural 
possibility that comes to mind is the Shannon entropy tr(pln p). However, this 
measure is sometimes inconvenient due to the appearance of the logarithm. 
Therefore, it is common to use instead the "hnear entropy" [23, 24, 25] 

S^l-tr{p') = ^{l-s-s), (18) 
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which is a good approximation for the Shannon entropy for almost pure states 
and is closely related to the so-called purity tr(p^). The history of this kind of 
measure is long, and can be traced back to at least 1939 [23]. The linear entropy 
possesses a convenient physical interpretation: Suppose that p describes an 
incoherent mixture of n orthogonal states with equal probability, then tr(p^) 
reads 1/n. For a single qubit the maxinmm linear entropy is therefore 1/2. It 
is zero if and only if p describes a pure state. 

The decoherence rate is well described by the entropy production 

S ^ -2ir{pp) ^ -s ■ k (19) 

which follows directly from the master equation, respectively from the equation 
of motion for the Bloch vector. Obviously, the coherent dynamics given by the 
first commutator in the master equation (12) does not increase the entropy. 

In the context of quantum computing, we are mainly interested in the entropy 
production for (almost) pure states, i.e. the initial entropy production which 
is determined by the eigenvalues (17) of the matrix M . Since these eigenvalues 
are invariant under unitary transformations, they represent a global measure 
for the infiuence of the dissipation. In particular, they are independent of the 
choice of the basis. In the worst case, the decoherence is determined by the 
eigenvalue with the largest real part. 

For the driven system considered below, the matrix becomes time-dependent 
and, consequently, the eigenvalues (17) must be replaced by the eigenvalues of 
a Floquet equation. Their computation, however, might be quite cumbersome. 
Therefore, a more convenient measure which is also applicable to driven sys- 
tems is the sum 7 of all eigenvalues which can be computed without explicit 
diagonalization by the trace of the matrix M, 

7 = trM = 2r. (20) 

Since 7 is real and larger than all real parts of the eigenvalues (which are 
positive), it represents an upper bound for the decoherence rate. Moreover, 

for the present case, the largest real part of an eigenvalue is at least 7/3 and, 
thus, we can conclude that 7 gives the correct order of magnitude for the 
decoherence rate of the most fragile initial state. 

A further related, but more probabilistic measure, that characterizes the loss of 
coherence, is the average S over all pure initial states of the entropy production 
and reads 

Tav = (^(0))av = |. (21) 
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3 Dissipation and driving 



To study the controllability of the decoherence, we generahze the master equa- 
tion to the case where a time-dependent Hamiltonian 

Hoit) = H^fit) (22) 

acts on the qubit. We assume that time-dependence f{t) is periodic with the 
driving period T = 27r/Jl and has zero average. Then, according to the Floquet 
theorem [26] , the corresponding propagator can always be written in the form 

U{t, t') = Up{t, t')UF{t - t') (23) 

where Up(t,t') = Up(t + T,t') obeys the time-periodicity of the driving field. 
The Floquet propagator Up depends only on the time difference and contains 
non-adiabatic phases which emerge during the propagation over one driving 
period [27]. Since Up is time-periodic, the long-time dynamics of the driven 
qubit is entirely determined by the one-period propagator UpiT). 

Inserting the propagator (23) into the master equation (12) results after the 
r-intcgration in a dissipative kernel that still depends periodically on the final 
time t. However, since we consider the case of fast driving, it is possible to 
separate timc-scalcs and replace the rapidly oscillating operators in the master 
equation by their time-averages. In doing so, we arrive at [22, 17] 

i 1 TTcy /\ 

P = -^[H^b + HD{t),p\ - -[a^, [Q,p\]+i^[a^,{ay,p}], (24) 

where the difference with respect to the master equation (12) for the static 
case comes from replacing in the second commutator Fa^ by the operator 

Q = -/ dt dTS{T)Ul{T)Ul{t-T,t)a,Up{t-r,t)UF{r). (25) 
1 Jo Jo 

The master equation (24) for the driven, dissipative qubit reflects the close 
resemblance to the master equation (12) valid in the static case. The last term 
in (24) is not modified by the driving. The reason for this is that the driv- 
ing enters linearly in (10) and, consequently, vanishes in the average over the 
driving period. However, besides the explicit presence of the time-dependent 
Hamiltonian in the coherent contribution, also the dissipative part has ac- 
quired a change: The coupling operator is now replaced by an operator 
Q that depends on the qubit propagator in the presence of the driving. As a 
consequence, we expect in the case of strong driving fields not only a modifi- 
cation of the coherent dynamics, but also of the coherence properties. In the 
following sections, we will investigate two typical types of driving: as a result, 
we find that they can alter the coherence times rather significantly. 
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4 Coherent destruction of tunneling 



As a first example for the significant influence of the driving field, we investi- 
gate the qubit under the influence of the Hamiltonian 

Hoit) = Aa:cCos{nt) (26) 

which couples like the bath to the qubit by via the operator ax and, thus, 
commutes with the qubit-bath coupling. Such a time-dependent field causes 
already interesting effects for the coherent qubit dynamics that we will briefly 
review before discussing decoherence. 

To sketch the coherent dynamics of the driven qubit, we first transform the 
Hamiltonian H^y, + Hoit) by the unitary operator 

Uo{t) = exp (-i^ sin(Oi)) . (27) 

This transformation defines the interaction picture with respect to Hoit) and 
results in the hkewise T-periodic Hamiltonian 

^qb(*) = ^ [(^z cos sin(Qt)) + ay sin sin(Qi)) | . (28) 

The corresponding Schrodinger equation cannot be integrated exactly since 
Hqh{t) does not commute with itself at different times and, thus, time-ordering 
has to be taken into account. We restrict ourselves to an approximate solution 
and neglect corrections of the order A^. Within this approximation, the propa- 
gator is simply given by the exponential of the integral of the time-dependent 
Hamiltonian. This consists of two parts: First, there is the time-average of 
Hqb{t) which determines the leading contribution to the Floquet propagator 
Upit). The second part comprises terms that oscillate with the driving pe- 
riod. For a high-frequency driving with Q ^ A, this latter contribution can 
be neglected since the periodic part of the propagator is dominated by the 
contribution Uq coming directly from the driving. This approach is equivalent 
to a perturbational computation of Floquet states [28, 29] as has been shown 
in [30]. 

Finally, the qubit propagator in the interaction picture defined by (27), is 
determined by the time average of the Hamiltonian (28) which reads 

^qb = ^CT,. (29) 

This static approximation to the driven qubit Hamiltonian is of the same 
form as the original Hamiltonian (2), but with the tunneling matrix element 
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renormalized according to 



A ^ AefF = Jo{2A/hn)A. (30) 

Here, Jq is the zeroth order Bessel function of the first kind. Thus, within 
this high-frequency approximation, the entire propagator for the qubit in the 
Schrodinger picture reads 

U (t, 0) = exp (^-i^ sm(nt)a^^ exp (^-i^j^a^^ . (31) 



Of particular interest are now driving parameters for which 2 A/ HQ corre- 
sponds to a zero of the Bessel function Jq, i.e. for which Aeg vanishes. Then 
the one-period propagator Uf{T), becomes the identity. Or in other words: the 
long-time dynamics is suppressed. The dynamics within the driving period re- 
quires a closer look at the periodic propagator Up: For an initial preparation 
in an eigenstate of it contributes only a global phase, such that the dy- 
namics as a whole is suppressed also within the driving period. This effect of 
suppressing the time-evolution by the purely coherent influence of an exter- 
nal field has been investigated first in the context of driven tunneling [13, 14] 
and is named "coherent destruction of tunneling" (CDT) ^ . However, for a 
preparation other than an eigenstate of a^-, the periodic propagator Up will 
still cause a non-trivial dynamics within the driving period. This fact has been 
found numerically in Ref. [10]. 

There is also a significant infiucncc of the driving (26) on the dissipativc dy- 
namics of a particle with an initial preparation in an eigenstate of (Tx'- Besides 
slowing down the coherent time evolution, also the dissipative transitions are 
decelerated, i.e. the rate for the dissipative transitions between the wells of a 
bistable system becomes lower. In the following, we investigate dccohcrence 
and dissipation under such a CDT driving without the restriction to a specific 
preparation. 

In order to evaluate the coefficients of the master equation (24) valid in the 
driven case, we have to compute the operator Q. This requires an explicit 
expression for the propagator U {t, t') starting at time t' and, thus, we have to 
consider initial phases. Repeating the calculation from the beginning of this 
section yields 

U{t — T,t) — exp ^— i^^[sin fl(t — t) — sin flt]ax^ exp ^— i — , (32) 

where the first factor is the time-periodic part Up{t — T,t) while the second 
factor determines the long-time dynamics. Since Up commutes with ax, only 

^ Although we consider here also parameters for which the dynamics is not entirely 
suppressed, we will refer to a driving of the form (26) as "CDT driving" . 
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the Floquet propagator Up is relevant for the operator (25). After performing 
the time and the frequency integration we find 



QCDT = ^«5(Aeff)(7a;. (33) 

Thus, the master equation is of the same form as in the undriven case, Eq. (12), 
but with r replaced by 

r.„.4.(A.,^r.,2./..)£^5~). (34) 

Since the spectral density S{u!) increases monotonically with the frequency 
cu and, moreover, the Bessel function Jo{x) < 1 for x > 0, dissipation and 
decoherence become smaller due to the driving. 

Again, we employ the trace of the matrix M appearing in the equation of 
motion for the Bloch vector as a criterion for the global decoherence strength 
and find 

7CDT = '5(Acfr) (35) 

which represents within the high-frequency approximation an upper hmit for 
the decoherence rates. For high temperatures, phA^s <^ 1, the coth in the 
power spectrum (14) can be replaced by the inverse of its argument. This 
defines the classical limit where the decoherence rate (35) becomes Ai^aksT /h, 
i.e. independent of the effective tunnel splitting Acfr. In the low temperature 
limit /3^Aeff ^ 1, the coth becomes unity and, thus, 

7CDT~7'/o(2A//ifi), (36) 

which means that decoherence becomes smaller by a factor Jo(2A//if2). This 
reduction of decoherence is brought about by the fact that the driving (26) de- 
celerates the long time dynamics of the qubit. Thereby, the frequencies which 
are relevant for the decoherence are shifted to a range where the spectral 
density of the bath is lower. Consequently, the infiuence of the bath is dimin- 
ished. Since 7cdt represents an upper bound for the decoherence rate, this 
demonstrates that coherence stabilization via CDT is a phenomenon that is 
independent of the initial state. 

It should be noted, however, that the coherent dynamics is slowed down by 
the same factor as the decoherence, cf. Eqs. (30) and (36). Thus, if for a 

specific application, the figure of merit is the number of coherent oscillations, 
the present coherence stabilization scheme may not prove useful. For resonant 
driving, this unfavorable situation can change for certain preparations, see 
Fig. 2 in Ref. [10]. 
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5 Dynamical decoupling 



A recently proposed mechanism for coherence stabihzation is the so-called dy- 
namical decoupling [18]. This scheme employs sequences of 7r-pulses which flip 
the sign of the operator that couples the qubit to the bath operators. The ba- 
sic idea originates from the suppression of spin diffusion in nuclear magnetic 
resonance experiments [31,32] and has become an established technique in 
that area [33]. In the present case where the bath couples to the operator 
[cf. Eq. (4)], such a transformation is e.g. induced by the Hamiltonian hcuRaz 
for a pulse duration tt/ujr. Since the corresponding propagator is a function 
of the qubit Hamiltonian, the coherent dynamics is not altered. Besides the 
prospective benefits of such a control scheme, there is also a number of possible 
drawbacks that the application of 7r-pulses might cause: For a driven system, 
there is always the possibility of unwanted off-resonant transitions [34] , espe- 
cially in the case of sharp pulses. A more practical limitation is the fact that 
only noise with frequencies below the pulse repetition rate can be eliminated 
in such a way. 

These disadvantages can be overcome partially by applying a continuous wave 
version of the dynamical decoupling scheme, i.e. a driving of the form 

Hd = Aa, cos{Qt) (37) 

for which the available frequency range is larger. This constitutes our second 
example of a driving field infiuencing considerably the qubit decoherence. In 

contrast to the driving Hamiltonian employed in the previous section, the 
present one commutes with the static qubit Hamiltonian and, consequently, 
the propagator for the driven qubit can be computed exactly to read 

U{t, t') = exp (-i^[sm{nt) - sm{nt')]a^\ exp (-iAcx^t - t')) . (38) 

Again, we have written the propagator in the form (23) which is suitable for 
simplifying the master equation (24). Inserting this into the expression (25) 
results for A ^ Q in the effective couphng operator 



Qdd = ^cj.. (^J^{2A/m)S{A) + '^p^Jn C^A/m) S{nQyj . (39) 



In order to derive this expression, we have decomposed the exponentials of the 
trigonometric functions into a Fourier series using the identity exp[ia; sin(fit)] = 
J2k Jkix) exp{ikQt) , where Jk is the kth order Bessel function of the first 
kind [35]. Like in the previous cases, the effective coupling operator Qdd is 
proportional to and, thus, the master equation is again of the form (12). 
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Fig. 1. Coherence stabilization factor r/ as a function of the driving frequency for 
dynamical decoupling for different temperatures. The cutoff frequency is lJc = 500A, 
2A/hn = 2.4, and the dissipation strength is a = 0.01. The dotted horizontal lines 
mark the values above which the coherence is improved on average (r/ > 1/4) and 
for an arbitrary initial condition (r/ > 1), respectively. 

The only change for the dissipative dynamics is the replacement of F by 



The decoherence rate in this case depends on the spectral density of the bath 
at multiples of the driving frequency which may be larger than the cutoff 
frequency Uc- The vr-pulses applied in the original version version [18] of dy- 
namical decoupling, correspond for a continuous driving to a field amplitude 
such that 2A/hQ equals the first zero of the Bcsscl function Jo. i.e. assumes 
a value 2.404825.... Then only the sum in Eq. (40) contributes to the de- 
coherence rate Fdd. If now the driving frequency is larger than the cutoff of 
the spectral density, fl > cUc, the decoherence rate is considerably reduced: 
For low temperatures, 1//3 <S ^A, the hyperbolic tangent in the decoherence 
rate (40) become unity and each contribution is weighted by a possibly large 
factor nQ/A. In the high-temperature limit 1/(3 <^ HQ, we use tanha; ^ x and 
find that the dependence of the prefactor on nfl cancels. This means that the 
dynamical decoupling scheme is especially useful for high temperatures. The 
physical reason for this is that the driving shifts the qubit dynamics towards 
high frequencies where the thermal occupation of the bath modes is negligible. 

To assess the coherence stabilization originating from dynamical decoupling, 
we define for the coherence stabilization the rather conservative measure 



F/2 



(41) 



7dd 
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i.e. the ratio between the lowest decoherence rate without driving and the 
largest decoherence rate with driving, cf. equations (17) and (20). For para- 
meters such that ?7 > 1, it is granted that the driving stabilizes the coherence 
independent of the initial state. If 77 > 1/4, we still can conclude that the 
coherence is improved in the average over all possible initial states. 

Figure 1 compares the coherence stabilization 77 as a function of the driving 
frequency for 2A/hQ = 2.4, i.e. close to a zero of the Bessel function Jq. It 
reveals that for driving frequencies well below the cutoff, the driving rather 
spoils the coherence. This improves with increasing driving frequency and, 
finally, for a high-frequency driving r) becomes much larger than unity corre- 
sponding to a significant coherence stabilization. The data demonstrate the 
usefulness of dynamical decoupling for high temperatures discussed above. 



6 Conclusions 



We have investigated the coherence properties of a qubit weakly coupled to a 
harmonic oscillator bath under the influence of an external driving fleld. Two 
types of driving have been taken into account: One that can cause coherent 
destruction of tunneling and another that corresponds to 7r-pulses that invert 
the sign of the qubit-bath coupling operator. To estimate the decoherence, we 
have derived an upper bound for the decoherence rate. As a result, we have 
found that both types of driving can enhance the coherence properties of the 
qubit significantly, independent of the preparation. 

In the case of CDT we have found that the main effect comes from the fact that 
the driving shifts the coherent long-time dynamics of the qubit towards lower 
frequencies where the spectral density of an ohmic bath is lower and, thus, the 
effective coupling is weaker. Consequently, decoherence is generally reduced. 
This is most pronounced at low temperatures. For high temperatures, however, 
the lower spectral density is counterbalanced by an increasing thermal noise, 
such that decoherence becomes independent of the driving. 

For the continuous-wave dynamical decoupling scheme, we have found that 
a low-frequency driving is counterproductive. However, once the frequency 
exceeds the bath cutoff, the coherence properties recover and are finally sig- 
nificantly improved, especially at high temperatures. Since such a dynamical 
decoupling by a harmonic driving allows higher driving frequencies than the 
pulsed version, this form of coherence stabilization bear interesting perspec- 
tives for applications. 
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